Thermodynamic approach to the dewetting instabiUty in ultrathin films 
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The fluid dynamics of the classical dewetting instability in ultrathin films is a non-linear process. However, 
the physical manifestation of the instability in terms of characteristic length and time scales can be described 
. by a linearized form of the initial conditions of the film's dynamics. Alternately, the thermodynamic approach 

7—i ' based on equating the rate of free energy decrease to the viscous dissipation [de Gennes, C. R. Acad. Paris. 

I v298, 1984] can give similar information. Here we have evaluated dewetting in the presence of thermocapillary 

OA ■ forces arising from a film-thickness (h) dependent temperature. Such a situation can be found during pulsed 

f— I ' laser melting of ultrathin metal films where nanoscale effects lead to a local h-dependent temperature. The 

I thermodynamic approach provides an analytical description of this thermocapillary dewetting. The results of 

i—i , this approach agree with those from linear theory and experimental observations provided the minimum value 

^ ' of viscous dissipation is equated to the rate of free energy decrease. The flow boundary condition that produces 

' this minimum viscous dissipation is when the film-substrate tangential stress is zero. The physical implication 

of this finding is that the spontaneous dewetting instability follows the path of minimum rate of energy loss. 

•1— I . 

o . 

C/3 , I. INTRODUCTION 

' Investigations of thin film morphology evolution and control is of fundamental and technological interest. In particular, 
^ ' spontaneous self-organizing processes \ Y] that lead to nanostructure formation in a reliable way have attracted tremendous 
^ attention. The resulting nanostructures can have novel behavior as well as be utilized in a wide variety of technologies, such 
as, energy harvesting |[2-01, biomedicineiH |^, spintronics |7], photonics |8] and magnetism IHITol]. One example of a self- 
H organizing process is the spontaneous dewetting of a continuous liquid film from a surface. A scientific understanding of 
dewetting has implications to many industrial applications, including in the deicing of airplane wings with non-wettable surfaces, 
in preventing hydroplaning of automobiles on wet roads due to thin continuous layer of water and in designing chemicals to 
prevent the break-up of the lachrymal film that protects the cornea of the eye. Another growing application of dewetting is in 
Q the fabrication of nanoscale structures in a robust, controllable and cost-effective manner The extensive studies of polymer thin 
films and growing number of investigations of dewetting in metalUc thin films is indicative of this technological interest and also 
to the need for a deeper understanding of the phenomenon lfTll - l20ll . 

The classical dewetting instability in thin films can be interpreted as a competition between two energy terms. For the case 
QQ of a large number of polymer or metallic films studied, these two energies correspond to the surface tension and the attractive 
intermolecular dispersion force between the film-substrate and film-vacuum interfaces mediated by the film material. As shown 
OO first by Vrij ll2lll22ll . the instability can be described from an energetic viewpoint by evaluating the thermodynamic free energy 
^""^ change of the system under perturbations to the film height. The prediction from such an energy analysis is that for certain 
. perturbation wave vectors, the film enters an unstable state and thus, can spontaneously dewet. As a result, studies of dewetting 
' have focused largely on the fluid dynamics of the film, through which it is possible to obtain the relationship between the rate 
0^ of growth or decay of surface perturbations to their wave vector, i.e. the dispersion relation. However, the fluid dynamics for 
even the simplest dewetting scenario, such as the example above, is a highly nonlinear process, and, while addressable by many 
J> [ numerical techniques ll23l - l26ll . is often evaluated through a linear analysis in order to achieve physical insights into dewetting. 
An alternate approach to quantitatively evaluate dewetting is thermodynamics. Fluid flow pathways can be analyzed through 
thermodynamic considerations in which the conversion^ useful internal energy to external energy loss via heat, such as by 
viscous dissipation, is used to quantify the behavior jZTlbsil . 

In this work, we show that such an approach can provide meaningful insight into the nature of fluid flow as well as the energy 
pathway for dewetting instabilities. Specifically, we have applied the thermodynamic formulation to the case of dewetting in 
which film thickness dependent Marangoni or thermocapillary forces are also present. Such a situation has been observed in the 
melting of thin metallic films by nanosecond pulsed lasers ll29il30ll . In our thermodynamic analysis the rate of thermodynamic 
free energy decrease due to film thickness fluctuations is balanced with the rate of energy loss due to viscous flow, i.e. viscous 
dissipation. This leads to a analytical description of the dewetting process without explicit need to solve the the height evolution 
dynamical equation from the Navier-Stokes (NS) equation. The thermodynamic and linear approach show identical results 
for classical dewetting. For thermocapillary dewetting, the two approaches agree only if the minimum viscous dissipation is 
evaluated. This minimum was found to occur for a particular flow boundary condition, which related the pressure gradient with 
the thermocapillary forces, and resulted in zero tangential stress at the film-substrate interface. Besides this physical insight into 
the fluid flow, the thermodynamic analysis also showed that the dewetting pathway is one in which the rate of energy loss is 
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minimized. 



11. THEORY 



For completeness, we first begin by summarizing the derivation of the thin film fluid velocity for a one-dimensional (ID) 
incompressible fluid from the NS equation within the lubrication approximation. A complete analysis is provided in ref. IBlll . 
In this approximation the average or unperturbed thickness /i„ of the film is much smaller than the in-plane dimension (x), as a 
result of which, the only velocity change of importance occurs along the thickness or z-direction. Furthermore, because of the 
small thickness of the film, inertial effects can be neglected and so the flow is dominated by the viscous effects. Using the above 
approximations, the NS equation for the steady-state condition (p ^ = 0) in the x-direction is given by: 

where v is the .ic-component of the liquid velocity, VP = ^ is the pressure gradient in the direction of flow x, and r\ is the 
dynamic viscosity. By integrating the velocity v as a function of height z we get: 

v^yp—+az + b (2) 
2tj 

Typical boundary conditions used to analyze the classical dewetting instability are the no-slip condition at the film-substrate, 
so v{h = 0) = 0, and a stress-free boundary condition at the top film surface. The no-slip condition results in b ~Q. At 
the top surface, we introduce the h-dependent Marangoni effect by equating the shear stress to the surface tension gradient: 

^5^l''o ~ ^'^inrix ~ y^'lx ~ iYh^')ho^ where, Y{h) is h-dependent surface tension of the film-vapor interface, | 7/, | = | ^ | is the 
magnitude of the height coefficient of surface tension, and =h' is the thickness/height gradient along the flow direction, with 
all quantities evaluated at the average film thickness ho. From this, the velocity and velocity gradient in the z-direction can be 
respectively expressed as: 

VP 2 VP/io - 7hh' 

V = —z z (3) 

2tj tj 

and 

dv VP VPho^Yhh' 

-J- — 2 (4) 

dz T] rj 

Based on this, one can now easily evaluate the rate of energy loss due to viscous liquid flow, i.e. the viscous dissipation, per unit 
volume e occurring in the film. This quantity is given by ll32ll : 

^ ^dz ri ^ TJ ^ 77 

For the case of classical dewetting, i.e. in which Marangoni forces are absent (jh = 0), the viscous dissipation will be: 

ec^^^^iz~h„)^ (6) 
J? 

where the superscript denotes classical. 

Next, we can evaluate the rate of thermodynamic free energy change AF for fluctuations/perturbations to the initial height of 
the film. Since we are primarily concerned with the dewetting instability, we will use the classical approach proposed by Vrij 
where-in the film-vapor surface tension energy competes with the attractive dispersion energy. Film height perturbations will 
increase the top film surface area and so surface tension increases the overall thermodynamic free energy of the film. On the 
other hand, the long range attractive dispersion energy varies as n ~A /2nh'^, where A is the Hamaker coefficient with negative 
sign, leading to an overall decrease in thermodynamic free energy. As Vrij showed, it is the free energy decrease resulting from 
competition between these two energy terms that drives the dewetting instability. Here, we evaluate the rate of this free energy 
change by expressing the height perturbations as Fourier components of type: 

h{x,t)^ho + £e'"e-"'''' (7) 
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where the perturbation has an amplitude of e, a characteristic temporal decay rate a and a corresponding wave vector k. In this 
work, we explore the thermodynamic solution for a thermocapillary/Marangoni problem where the surface tension is dependent 
on the local height but independent of the x-position on the surface. Consequently, the surface tension change does not contribute 
to a change in free energy. This can be shown as follows. The standard procedure is to calculate the rate of total free energy 
change due to the height perturbation and evaluate it over one wavelength of the perturbation Ii2li1 . In this case, the resulting rate 
of energy change will be 



AFs 



single wavelength 



•1 .dh 



-n Ah^ 

2 



■ -YhAh)dx 



(8) 



One can immediately see that when the magnitude of the surface tension derivative with respect to height (7/, = dy/dh) is 
independent of the x-position along the film, the integral of the third term is Jg 5 ^£e"'*^^+'^'t/ji: = 0. Consequendy, to calculate 
the rate of change in the free energy, we can ignore the contribution from 7/,, and express it as the difference due to the initial 
and perturbed film thickness at any position x 112211 as: 



dt 



1 , dh 1 // . , 9 



2nK 



2 2{at-ikx) 



(9) 



The first term in the expression on the right hand side is the rate of increase of surface tension energy, and the second one is the 
rate of change in energy from the dispersive interaction. 



III. RESULTS 



A. Dispersion relation for classical dewetting 



Here we compare the characteristic dewetting length scales obtained from fluid dynamics versus the thermodynamic approach. 
The typical approach to obtain the dispersion relation between the rate o and wave vector k has been to describe the rate of 
change in film height based on the NS equation and mass conservation ll2ll l22ll. For the classical dewetting instability, i.e. 
without Marangoni effects, the resulting dynamics is described by the equation lll8l[33ll : 



dh f A 

3ii— = -V-(7h'v-\7^h- — ^h] (10) 

As is evident, this equation is non-linear in h and presents considerable challenges towards achieving an analytical description 
of dewetting that could provide simple but physically insightful information about the instability. Consequently, a prevalent 
approach is a solution afforded by linear stability analysis. Hence, after an expansion of the right hand side of Eq. [TOl 

= (-/t2£e^'-''-)f7fc2(/,3 + 6/i2ee<^'-*- + 9/i„eV('^'-*'=)+4£3^ — — — (11) 
dt \ 27t{ho + ee"' iKxy J 

and applying by keeping terms only linear in the perturbation amplitude £, as required by LSA, Eq. [TT]reduces to : 

3r]{eoe'"-'^') = {^khe""-'^') (rk^hl + (12) 



Consequently, the resulting dispersion relation is given by ||29|] : 

<^c-^-f»^ + 4) (.3) 

The characteristic (or classical) dewetting length scale A^'^'^ can be obtained from the dispersion relation by the maxima 
condition ^ = and leads to: 

^u^jjji^f^n^,. (.4, 
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where the superscript refers to LSA. 

On the other hand, the thermodynamic (TH) approach is based on equating the rate of free energy change (Eq. |9]l to the total 
viscous dissipation in the film. We can calculate the total viscous dissipation per unit area E for the liquid film by integrating 
over the film thickness as follows: 

Ec = jecdz = ^^^hl (15) 



The next, and important, step in evaluating this integral is to relate the pressure gradient to the film height through a volume 
conservation argument. Volume conservation requires that the rate of change of film height dh jdt be related to the flux of liquid 
flow J(x) as dh/dt — —V.J{x). To evaluate this we have used the thin film lubrication approximation in which flux can be 
written in terms of the height-averaged liquid velocity < v > as [31]: 

1 VP , 

J{x)=ho<v>=ho.{— vdz)^-—K (16) 
ho Jo 3t] 

from which we can express the volume conservation equation as: 

dh „ V2p 
dt 377 

By rearranging terms we get the desired relation for the pressure gradient as follows: 

VP^ [ V^Pdx - i (l!?^)ee^'-'^- (18) 
J k h^ 

where we have made use of Eq. |2l Using this expression in Eq. [15] the total dissipation is: 

£c = -4^cj2eV(^'-''^^') (19) 

Then, equating Eq. [T9]with the rate of decrease of free energy, Eq. |9] AF = Eq, and expressing the result in terms of a we get 
the dispersion relation from the thermodynamic approach: 



from which, the classical dewetting length scale can be expressed as: 

1/2 



TH , 



,r„IJ_^\ ,21, 

where the superscript and subscript refers to thermodynamic. As expected, the fluid dynamics described by LSA and the TH 
approach give identical results for the classical dewetting instability jZTllsill . On the other hand, as we show next, dewetting 
with Marangoni (or Thermocapillary) forces requires a more stringent evaluation of the viscous dissipation in order to provide 
results comparable to LSA. 



B. Dispersion relation for thermocapillary dewetting from thermodynamics 

As in the previous section, we analyze the LSA and TH approaches for dewetting in the presence of thickness-dependent 
Marangoni effects. In order to relate our work to experimental observations, we describe LSA results for the case when ultrathin 
metal films on Si02 substrates are melted by nanosecond ultraviolet wavelength laser pulses liTii ITtIi . In this situation, there 
is a strong thickness-dependent reflection and absorption of light by the thin metal film which leads to a local h-dependent 
temperature of the liquid metal film ifTgll . In addition, as reported previously, the temperature gradient along the plane of the 
film, dT /dx, generated by this nanoscale heating effect, can have a positive or negative sign depending upon the initial film 
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thickness ho II29I1 . With this, the boundary condition describing the /z-dependent Marangoni effect can be rewritten in the form 
of a thermocapillary effect as follows: 

Where ^ is the temperature coefficient of surface tension, Tf, — dT /dh is the film height-dependent temperature and 
h' — dh/dx. Given that all metals have a negative value of Jt, we have expressed the boundary condition in a more conventional 
form using 77- = — | I- In this scenario, the resulting dynamical equation of the film height is given by: 

dh / , 9 A 3 \ 

3^7 ^ = - V • ( 7/1' V • V' - ^ V + 2 ^^rj (23) 

Applying steps similar to to the classical case, the right hand side of Eq. l23]is: 

3t]|J = - {ee"' {yk^ho + ee"" f [ho + Aee'"-'''''^ (24) 



2n{ho + ee'"-''^)^ ' " ' " '^2 2 

which, upon applying LSA, reduces to: 



and the resulting dispersion relation is ll29ll : 



3r] ^'^ ^iKhyi ho 



^Tc' - -^irk' + + ^^^) (26) 



From Eq. |26]the characteristic dewetting length scale in the presence of thermocapillary effects, A^, can be expressed as: 

y/^ 2 

V A + 3n\rT\hlT, 



where the superscript TC denotes thermocapillary. 

Next we evaluate the length scale using the TH approach based on evaluating the the total viscous dissipation per unit area E for 
the liquid film. First, the dissipation per unit volume expressed in terms of the thermocapillary boundary condition is: 

e7-c = J?(^) =—^\z-ho) -2VP [z-ho)^ (28) 

The total viscous dissipation per unit area E for the liquid film is now: 



Er, = Lc* = + V/.^^*; + <li^^*, (29) 

J 3r\ " rj T7 



Again, the next step in evaluating this integral is to relate the pressure gradient to the film height through volume conservation 
arguments, as done for the classical case. 

Ji.)^ho<.>^ho.il- f\dz)^-^-^hl-^-J^Ll^ (30) 

from which we can express the volume conservation equation as: 

dh V^f 3 \7r\T„h"hl 
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By rearranging terms and substituting the height perturbation, Eq. |7] into above we get: 

The pressure gradient along the x-direction can now be obtained by integrating Eq. [32]as follows: 

VP = [ V^Pdx = A + (33) 
J k hi 2ho 

Then, on substituting the above relation for pressure gradient into Eq. |29] we get the total viscous dissipation: 

p / pj? ^2^l\rrVh^ , 3(1 yt I nfk^h,, 



hlk"^ ho 477 

^^3r^\rT\n^^3i\rT\nn^^hJjro^^ ^ I r.)VUv(— ) (34) 
«o 2ho ri ri } 

Finally, equating Eq. [34] with the rate of decrease of free energy, Eq. |9]and rearranging the equation in terms of CT we get an 
analytical dispersion expression for thermocapillary dewetting as: 

It is important to enote that Eq. [35l which comes from the thermodynamic approach, differs from the LSA result, Eq. |26l One 
can note that there is substantial difference between the two solutions. The TH approach is quadratic in a as well as | ■j^ | T/, 
(Eq. [35] l. while, in the LSA case (Eq. |26] | it is linear in both quantities. Consequently, the TH approach does not directly 
lead to the linear dependence on | ■j^ | Ti, as evident from LSA (Eq. |26] |. This is especially important since, as noted earlier, 
the thermal gradients generated by pulsed laser heating can have positive or negative signs and hence lead to fundamentally 
different dewetting behaviors [29, 30]. Since the TH dispersion is a quadratic function of the thermal gradient, its behavior will 
be independent of the sign of the thermal gradient and so does not agree with LSA. As we show next, it is necessary to evaluate 
the characteristics of dissipation in order to get the correct behavior from TH. 



1. Dispersion using minimum viscous dissipation 



In the classical case, the total viscous dissipation is uniquely defined by the magnitude of the pressure gradient for any given 
film thickness, as evident from Eq. [15] On the other hand, the total dissipation for thermocapillary dewetting is not unique, 
and, in fact, varies with the magnitude of the thermal gradient for any given pressure gradient, as evident from Eq. [29] It is this 
behavior that is responsible for the above discrepancy between LSA and TH and can be resolved by evaluating the minimum 
viscous dissipation. 

The minimum viscous dissipation for the fluid being subjected to pressure gradients can be estimated from the differential 
condition de/dVP = 0. Using Eq. [28] this leads to the condition Vf (z — ho) =| 'j^ | 7/,/i'. The general solutions satisfying the 
above equality can be evaluated for various values of height z in relation to the thickness h^. First, the condition z = ho does not 
yield a unique relation between VP and | 7/, | h' and is therefore not a useful solution in the context of the dissipation. On the 
other hand, the choice of z = 0, yields the case of VP - . By utilizing Eq. [4] one can see that the physical interpretation 
of this condition is that the tangential stress at the film-substrate interface at z is zero. The resulting viscous dissipation for 
thermocapillary dewetting will now be (from Eq. l28l i: 

^TC = — - — (36) 

where the superscript m signifies a minimum. One can verify that this is a minima by noting that the second derivative d^Cy/ d^VP 
is positive. Therefore, the minimum dissipation per unit area of the film can be obtained as: 

E¥c = Je7cdz=^^^ (37) 

o 
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?^(/ + ^)±^\/(/ + ^)^-^^ (40) 



Using the form of Vf from Eq. [33] the relevant form of fi™^ is: 

Finally, by equating the rate of free energy change AF (Eq. |9|l and the minimum viscous dissipation rate E^c (Eq. [38T l. we obtain 
an analytical form of the dispersion relation as: 

This result is clearly different from Eq. [35]because now, the linear behavior with | 7r | 7/i is also present. Next, we evaluate this 
quadratic equation in o for various magnitudes of the thermal gradients and show that it is identical to the LSA results for large 
thermal gradients. By defining / = yk^ +A/2nhQ and ^ = (3 | ■j'r | Ti,)/ho, the roots of the dispersion relation are: 

1. Minimum dissipation approach for classical case (i.e. = 0) 

By substituting Ti, = g = in Eq. |40l the resulting relevant root is: 

As expected, this result is identical to the result for the classical dewetting instability. 

2. Minimum dissipation approach for strong thermal gradients 

In the case when the magnitudes of the thermal gradients are larger then the attractive dispersion forces, i.e. for example 
when 3 | 7^ 7), | > | A /27:hl \ , then we have | g | > | / | . This is the condition found in the experimental cases reported earlier 
ll29l[30ll . and we can obtain an approximate solution from Eq. |40]as follows: 

<..S-^(/ + .)±^v^=-|^(/ + ri (42) 

or 



From Eq. |43] the characteristic wavelength Aj^ obtained from the maxima in the dispersion given by da/dk = is: 



Comparing with the LSA result, Eq. |27], the only difference is in the factor of two multiplying the thermal gradient term. 
More importantly, the functional dependence on /!„, A, 7 and | 'j^ | 7/, remains the same, confirming that the minimum 
dissipation approach gives similar physical characteristics for the dewetting instability. An important benefit of utilizing 
the thermodynamic approach is also evident here. From the above dissipation analysis it is clear that there are multiple 
choices for the dewetting pathway in regards to the rate of dissipation. However the instability clearly picks the path 
which minimizes the rate of this dissipation, or in other words, minimizes the overall rate at which energy is lost in the 
dewetting process. We have also plotted the behavior of the dispersion relation, Eq. |43]for Co films on Si02 substrates 
with the various materials parameter values being: Hamaker coefficient A = —1.41 x 10^'** h J — 1-88 J/m^ and jj = 
—4.8 X lO^^y/m^ — K. Fig. Wia) plots the growth rate a versus wave number k for the dispersion relation Eq. |43] The 
data was evaluated for Co film of thickness 8 nm and various values of the gradient 7},, including T/, = 0, T/, > and 
Tf, < 0, as indicated on the plot. In Fig. [Tfb) ) we have plotted the characteristic length scale Aj^ for Co films on Si02 
as a function of varying thickness and various T/,. The decrease in length scale with increase in the magnitude of 7), < 
is evident. In Fig. [Tfc) the cut-off wavelengthA^^^ ^_,^y as a function of varying thickness and various 7), is plotted. The 
cut-off wavelength corresponds to the intersection of the growth rate with the wave number axis in Fig. [TJa). 
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3. Minimum dissipation approach for weak thermal gradients 

In the case when the magnitudes of the thermal gradients are smaller then the attractive dispersion forces, i.e. for example 
when 3 | JrTh \<\ A/lTihl |, then we have | g |<| / |. In this situation, we can approximate Eq.|40]as follows; 

where we have used the binomial approximation (1 — [g/f YY^'^ ^ (1 ^^^/2/^). Here we find again that the dispersion is 
considerably different from the LSA result of Eq. |26l 



IV. CONCLUSION 



We have theoretically evaluated the classical and thermocapillary dewetting instability in thin fluid films via a thermodynamic 
approach. In this, the rate of change of free energy is equated to the viscous dissipation in the thin film. The thermodynamic 
approach leads to an analytical expression for the dispersion without the need for explict solutionof the height evolution dynam- 
ical equation from NS within the lubrication approximation. We have compared results from this approach to existing results 
obtained by linearization of the fluid dynamics of the thin film. For the case of classical dewetting in the presence of surface 
tension and long range attractive forces, the thermodynamic approach predicts identical behavior to that from linear analysis. 
We have also evaluated dewetting in the presence of film-thickness dependent temperature variations. Such a situation can be 
found during dewetting of thin metallic films melted by a nanosecond ultraviolet pulsed laser. In this condition, a film thickness 
dependent reflection and absorption leads to thermocapillary forces along the plane of the film. In this scenario we found that 
the thermodynamic approach agrees with linear analysis provided the minimum viscous dissipation is evaluated. The fluid flow 
condition that gives minimum viscous dissipation is one where the film-substrate tangential stress is zero. In the context of 
dewetting in the presence of film thickness dependent thermocapillary forces, the thermodynamic approach clearly illustrates 
that the instability chooses a pathway which minimizes the rate of energy loss in the system. This results shows that the thermo- 
dynamic approach based on evaluating the rates of free energy change and energy loss is a simple but potentially powerful way 
to gain physically meaningful insight into such spontaneous pattern formation processes. 

RK acknowledges support by the National Science Foundation through CAREER grant NSF-DMI-0449258, grant NSF- 
CMMI-0855949, and grant NSF-DMR-0856707. 
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Figure captions 

• Figure [H (a) Plot of the growth rate versus wave number for the dispersion relation obtained from the dissipation 
approach for strong thermal gradients. The data was evaluated for Co film of thickness 8 nm and various values of the 
gradient Ti„ as indicated on the plot, (b) Plot of the characteristic length scale Aj^ for Co films on Si02 as a function of 
various T^. (c) Plot of the cut-off wavelengthAcutt-off for Co films on Si02 as a function of various Tf,. The cut-off wave 
number corresponds to the intersection of the rate with the wave number axis in Fig. (a). 
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